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A  simple  model  is  presented  which  yields  analytical  expressions  for 
the  heating  efficiency  of  beat-wave  excitation  in  a  plasma  with  a 
linear  density  profile.  The  effect  of  self-consistent  Landau  damping 
by  tail  electrons  is  included  without  recourse  to  WKB  approximations. 
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Presently  there  Is  renewed  Interest  in  beat-wave  excitation1-4  of 
Langmuir  waves  to  generate  fast  electrons,  two  possible  applications  being 
advanced  particle  accelerators5  and  modification  of  the  earth's  ionosphere. 
Although  the  description  of  this  process  in  uniform  plasmas  is  relatively 
well  understood,  the  presence  of  a  density  gradient  introduces  several 
complications  which  have  been  handled  in  a  WKB  sense  in  the  pioneer  work  by 
Rosenbluth  and  Liu2  and  in  computer  simulations  by  Cohen,  et  al4.  Because 
in  certain  experiments  the  role  of  density  gradients  may  not  be  negligible 
or  in  fact  may  be  useful  in  achieving  a  desired  result,  it  is  convenient  to 
obtain  analytical  expressions  which  explicitly  exhibit  the  dependence  of  the 
heating  efficiency  on  the  density  scale  length  L.  In  this  brief  note  we 
present  a  relatively  simple  model  which  yields  such  a  result  and  incorporates 
the  self-consistent  effect  of  Landau  damping  by  fast  electrons  without 
recourse  to  WKB  approximations.  It  is  hoped  that  this  result  may  be  useful  in 
estimating  the  performance  of  future  experiments  and  that  the  model  may 
provide  a  basis  for  more  elaborate  calculations  of  analogous  phenomena  in 
nonuniforra  plasmas. 

The  one-dimensional  model  consists  of  separating  the  roles  of  warm 
background  electrons  having  a  spatially  varying  density  n0(z),  and  fast  tail 
electrons  whose  density  n^  is  considered  to  be  uniform  (because  of  their 
longer  collision  mean  free  path)  over  the  region  where  the  relevant 
wave-particle  interactions  occur.  The  tail  density  is  assumed  small  (nt/n0  << 
l)  so  that  the  principal  contribution  of  the  fast  particles  enters  through 
Landau  damping,  while  the  warm  background  electrons  determine  the  propagation 
features  ;f  rhe  driven  Langmuir  wave.  The  4s^kgtou.iu  c lecL >itc  LLcated  cs 
a  warm  fluid  with  thermal  velocity  v  and  the  tail  electrons  are  described 
kinetlcally  assuming  an  unperturbed  tail  distribution 
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:ot  =  T"2Vvt 2 )  1  / 2  exp[-  v  /2vt  ] 


Beat  excitation  is  envisioned  to  arise  from  the  beat  ponderomotive  force 
acting  on  the  background  electrons  and  caused  by  two  transparent  electro¬ 
magnetic  waves  having  frequencies  and  wave  numbers  cjj  ,  kj  ;  j  =  1,2.  These 
waves  propagate  along  the  nonuniform  direction  z  and  can  individually  move  in 
the  direction  oi  decreasing  density  (kj  >  o)  or  towards  increasing  density 
(kj  <  0).  The  self-consistent  beat  excited  longitudinal  electric  field  has  a 
harmonic  time  dependence  E(z)  exp(-imt)  with  m  =  w  -  m  ,  and  is  determined 
from  Poisson's  equation 

"57  fE(z)  ”  Eo(z)  J  =  -4jre  [n  (z)  +  nt(z)j  ,  (2) 


in  which  the  effective  pump  source  E0(z)  for  beat  excitation  is 


E°  -  -  E«>  E< 


where  E0j  is  the  complex  amplitude  of  the  jth  electromagnetic  wave,  e  and  m 
are  the  charge  and  mass  of  an  electron,  and  c  is  the  speed  of  light.  In  Eq . 
(3)  the  frequencies  of  the  electromagnetic  waves  are  assumed  high  enough  that 
their  phase  velocities  are  essentially  the  speed  of  light.  The  perturbed 
charge  densities  (oscillating  at  frequency  u>)  for  the  background  n0  and  tail 
nt  are  calculated  from  the  fluid  and  Viasov  equations,  respectively. 

Introducing  the  Fourier  transform  of  the  driven  field 


E(z)  -  ~  /  dk  E(k)eikz  , 


and  assuming  a  linear  density  profile  with  scale  length  L,  i.e.. 


1  -  u)p2(z)/u>2  =  z/L,  where  u)p(z)  is  the  local  electron  plasma  frequency, 
transforms  Eq.  (2)  into  a  differential  equation  in  k-space 

rik  E(k)  +  [ifkS(k)  -  kp-]E(k)  =  2irHo«(k  -  kb)  .  '.5) 

with  kj,  =  k2  -  k;^  the  beat  wave  number,  kp  =  w/v  ,  and 

g(k)  =  /7  (  (  ~  }  kD  sgn(k)expf-  -~r~t >  (6> 

where  sgn(k)  =  1  for  k  >  0,  -1  for  k  <  0. 

Equation  (5)  can  be  integrated  to  yield  the  spectrum 

E(k)  =  -  2niEoL0(k  -  kb)exp[F(k)  -  F(kb)]  ,  (7) 

F(k)  =  -i  -  g(k)L  ,  (8) 

kD 

where  9(k)  is  the  Heaviside  step  -  function. 

The  sign  of  kb  determines  the  direction  in  which  the  directly  driven 

Langmuir  wave  propagates.  For  kb  >  o,  propagation  occurs  in  the  direction  of 

decreasing  density  and  results  in  a  single  wave  whose  phase  velocity  decreases 

_  _  -1/2 

as  z  increases  (in  the  WKB  sense  as  /3  v  (z/L)  ).  For  kb  <  0  the  spatial 

pattern6  consists  of  two  waves.  The  directly  driven  wave  propagates  towards 
increasing  density  and  upon  reaching  the  cut-off  point  (z  =  0  where  ojp  =  w) 
exhibits  partial  (because  of  Landau  damping)  reflection  and  thus  generates 
a  second  wave  which  behaves  analogously  to  the  one  directly  driven  for  kb  >  0. 
These  features  imply  that  for  kb  >  0,  tail  electrons  with  velocity  v  >  0  are 


■  ■  V  *.  - 
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accelerated  and  this  results  In  an  increase  of  the  heat  flux  SQ+  as  z  +  ■*>. 
However,  for  k^  <  0,  tall  modification  occurs  for  both  v  <  0  and  v  >  0 
particles,  and  causes  heat  flux  modifications  in  both  directions.  The 
self-consistent  modifications  in  the  time-averaged  tail  distribution  function, 
accurate  to  second-order  in  the  amplitude  of  the  longitudinal  electric  field, 


<«t  (z  .  i  v»  -  (  V>2  V17<  lE~  <“  -7  >  I2  H  4  t  • 


where  the  label  ±  refers  to  particles  moving  in  the  positive  or  negative  z 
direction. 

In  determining  the  efficiency  of  energy  transfer  to  the  plasma  the 
quantities  of  interest  are  the  asymptotic  heat  fluxes 


±co  mv3 


r  ,,  .  in  v 

5Q+  =  ±  f  dv  — T 
o  2 


<  6f  > 


which  using  Eq .  (9)  and  integrating  by  parts  can  be  put  in  the  form 


=  t  £  />  |E(k 


3v  ot 


From  Eq . ( 6 )  it  is  seen  that  g(k  =  w/v)  is  proportional  to  fot(v)  and  hence 
|  E  (k  =  w/v) | 2  depends  exponentially  on  fot(v).  This  implies  that  the 
integral  in  Eq .(11)  can  be  evaluated  analytically  without  approximations. 

Defining  the  heating  efficiency  n+  as  the  ratio  of  the  enhanced  tail  heat 


flux  to  the  Poyntlng  flux  Pj  =  c  |  E  j  [  2/8ti  of  one  (say  j  =  1)  of  the 


electromagnetic  waves,  it  is  found  that  for  direct  excitation  with  k^  >  0 
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5Q+  -  w5  L  P2 

n+  “  Pj  _  io |  to|  (c>(  n0 


where  a  =  (tt/2)i/2  (nt/n0)  (v/vt) 
density  and  scale  length,  and  £  = 
density  nG  is  evaluated  at  z  =  0. 


—  3)  I  1  -  exp[-2a( 1 
me 0  1 

(k  L)  measures  the 
D 

w/(/2  k  vt).  In 
D  c 


-  e~?2)]l  ,  U2) 

combined  effect  of  tail 
Eq .  (12)  the  background 


In  the  limit  of  large  a(i.e.,  long  scale  length) 
(o  ?2 

n+  a  2ir  (^)S  (k°JL)  (  n^tnc 2 ^  "  710 


(13) 


where  k0j  =  mj/c  since  oij  »  u)2  »  to.  It  should  be  noted  that  in  this  limit 
the  heating  efficiency  is  independent  of  tail  parameters.  In  the  opposite 
limit,  i.e.,  a  <<  1,  corresponding  to  small  tail  densities  and/or  short  scale 
lengths , 

n+  =>  n0a(  1  -  e  ~Z2)  ,  (14) 


showing  a  scaling  proportional  to  nt  and  L2. 


For  excitation  towards  increasing  density,  i.e.,  k^  <  0,  the 
contributions  from  the  two  heat  fluxes  at  z  =  ±  <»  must  be  included  in 
evaluating  the  efficiency,  which  yields 

_  6Q+  -  5Q- 

n-  =  - ^ 

=  n0  { 1  “  exp[-2a(l  +  e~^ 2 ) ] }  ,  (15) 

and  for  large  a  reduces  to  n-  »  Cq,  as  in  the  kb  >  0  case.  In  the 
small  a  limit  Eq.(15)  becomes 


k.  <•.  »"■  ^  «*-  " 


n0ct(  1  +  e“s  ) 


For  transparent  electromagnetic  waves  propagating  in  the  same  direction 


>>  1.  In  this 


(i.e.,  k_  ^ .  k_2  >  0)  it  is  likeLy  that  in  many  applications  £ 
special  case  it  is  found  that  <5Q_  +  0,  and  thus 


n+  =  n_  =.  n0 (  l  -  exp(-2n)) 


The  physics  behind  Eq .(17)  is  that  particle  acceleration  results  from  the 
transit-time  of  v  >  0  particles  through  the  spatially  varying  envelope  (and 
not  the  phase)  of  the  beat  excited  wave. 


In  summary,  a  simple  model  has  been  developed  that  yields  a  compact 
analytic  expression 


n±  =  no  { i  ~  exp(-2a(I  +  e  =2)J 


for  the  heating  efficiency  of  beat-wave  excitation  in  the  direction  of 
decreasing  (  +  )  and  increasing  (-)  density  in  a  nonuniform  plasma.  This  result 
includes  nonuniform  wave  propagation  and  self-consistent  Landau  damping 
without  recourse  to  WKB  approximations. 
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